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We introduce an extension of a variationally optimized perturbation method, by combining it 
with renormalization group properties in a straightforward (perturbative) form. This leads to a 
very transparent and efficient procedure, with a clear improvement of the non-perturbative results 
with respect to previous similar variational approaches. This is illustrated here by deriving opti- 
mized results for the mass gap of the 0{2N) Gross-Neveu model, compared with the exactly know 
results for arbitrary A''. At large A'", the exact result is reproduced already at the very first order 
of the modified perturbation using this procedure. For arbitrary values of N , using the original 
perturbative information only known at two-loop order, we obtain a controllable percent accuracy 
or less, for any A'^ value, as compared with the exactly known result for the mass gap from the 
thermodynamical Bethe Ansatz. The procedure is very general and can be extended straightfor- 
wardly to any renormalizable Lagrangian model, being systematically improvable provided that a 
knowledge of enough perturbative orders of the relevant quantities is available. 
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I. INTRODUCTION 

The variationally improved or optimized perturbation (OPT) is by now a rather well-used modification of standard 
perturbation theory (for a far from complete list of early references, see e.g. [1]). It is based on a reorganization 
, of the interacting Lagrangian such that it depends on an arbitrary (mass) parameter (so-called linear 5 expansion 
(LDE) in its simplest form), to be fixed by a definite optimization prescription, but it has many other variants [2-4]. 
In _D = 1 theories, such as the quantum mechanical anharmonic oscillator [5], described by a scalar (f)^ field theory, 

■ the LDE turns out to be very similar [6] to the "order-dependent mapping" (ODM) resummation method [3], being 
. equivalent at large orders to a rescaling of the adjustable oscillator mass with perturbative order, which can suppress 

the factorial large-order behaviour of ordinary perturbative coefficients. This appropriate rescaling of the adjustable 
^ mass gives a convergent series [6, 7] e.g. for the oscillator energy levels [5] and related quantities. For the oscillator 
field theory, no renormalization is needed, moreover the (ordinary) perturbative series is known to arbitrary high 

■ orders, and the known large order behaviour of the series is a crucial guide both numerically and analytically to 
. . , establish such convergence properties. In contrast, for most D > \ models, things become more involved due to 

. 5^ ' the necessary renormalization, but the OPT procedure can be made fully consistent [8] with the renormalization 
' program of ordinary perturbation theory, at least in a minimal subtraction scheme. Accordingly, the prescription 

■ is well-defined on renormalized Lagrangians with appropriate countcrtcrms. It is such that any physical quantity 
whose ordinary renormalized perturbative series is available can then be evaluated to order S'^ using well-defined 
modified Feynman rules. However in most D > 1 renormalizable models the perturbative series is only available for 
a few first orders, such that one can hardly study convergence properties of the OPT method. Nevertheless it gives 
at least a well-defined systematically improvable way to go beyond mean field approximation and has a wide range 
of applicability, also at finite temperature and density. To quote just a few examples of successful applications, for 
instance the occurence and precise location of a tricritical point and mixed liquid-gas phase within the GN model [9] 
in D = 2 -I- 1 dimensions [10] have been obtained, that was hinted at by Monte-Carlo simulations [11] but completely 
missed by the mean field approximation. OPT has also been used very recently [12] in the non- renormalizable 
_D = 4 Nambu-Jona-Lasino model [13] to study the phase transitions beyond mean field approximation in this 
simplified picture of low-energy QCD. In a different context, results for the shift of the critical temperature due to 
interactions in the Bose- Einstein condensate have been obtained [14-19], with some results [18, 19] in remarkable 
agreement with precise results from Monte-Carlo lattice simulations [20]. In this latter case the relevant field model 
is the three-dimensional super-renormalizable 0{N) model [21-23], which somewhat simplifies renormalization 
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issues, and the perturbative series is known to high (seventh) order, which also allows some definite conclusions on 
convergence properties[18, 19]. 

Now, in a more general renormalizable theory, once the renormalization procedure is well-defined, it is highly 
desirable to examine how the renormalization group (RG) invariance of the full theory is realized within this optimized 

perturbation. More practically, it is also of interest to examine how to incorporate eventually more information on 
higher orders from the RG within the calculated quantities. However in order to keep the RG resummation full 
information within the delta-expansion, it turns out that one has to resum to all orders in the latter expansion. 
This can be done [8, 24-26], at the same level of approximation than RG is treated {i.e. typically for the leading, 
next-to-leading, etc. . . logarithmic dependence on the mass parameter), but implying a rather involved formalism, 
rendering the practical optimization with respect to the perturbative mass parameter not much straightforward 
(and only numerical). Numerical estimates in reasonable agreement [8] with the exactly known mass gap of the 
0{2N) Gross-Neveu (GN) model have been obtained in this RG-resummed 5-expansion (at the few percent accuracy 
level). The method has also been applied to the QCD basic Lagrangian, and results have been obtained for relevant 
chiral symmetry-breaking quantities like the (constituant) quark mass, the quark chiral condensate and pion decay 
constant [25, 26], using only the (original) perturbative informations available at second order for those quantities. 
However, admittedly there is not much evidence for a systematic control of the convergence and error of the method 
from those results, even e.g. in the GN case where a precise quantitative comparison with known exact mass gap 
results for arbitrary N values can be addressed. For instance in the GN model, the discrepancy with the exact mass 
gap M{N) depended quite a lot on N, and also on the practical method used to optimize the mass parameter-'^. 
It was also shown [27] that at large perturbative orders the method gives a damping of the generic factorial growth of 
the large order perturbative coefficients (those due to the usual renormalons), but this modification is not sufficient 
to establish formal convergence of the new series, due to the extra complication brought about by logarithmic 
dependence in the mass in such renormalizable theories. A Borel resummation can improve the situation [27], 
but in any case such essentially qualitative large order results are of limited practical use to determine precise 
non-perturbative predictions for relevant physical quantities, for which in most cases only the very first few order 
perturbative coefficients are known exactly. 

Despite the lack of a rigorous convergence proof of the method for arbitrary renormalizable models, one may try to 
improve the basic method in order to obtain eventually more efficient prescriptions and better approximations. In this 
paper we reconsider the basic LDE-OPT construction, but augmented from the very beginning by the information on 
RG invariance properties of the physical quantities that are calculated. More precisely, we require the RG invariance 
to hold at the purely perturbative level of the 5-expansion and in a most straightforward manner i.e. consistently at 
the same order as the perturbative (non-RG) available information. This is in a sense a less ambitious program than 
constructing explicitly RG-resummed quantities, but as we will sec it has the advantage of giving a very transparent 
procedure with RG properties 'enrooted' at any step, and using nothing else than the available original purely pertur- 
bative information. This new procedure has immediately as consequence to link the mass and coupling parameters, 
in a way similar to the standard RG behaviour. But in addition, since optimization also fixes the perturbative mass 
parameter, it implies that both the mass and the coupling are fixed, and there are no free parameters, for any given 
values of the optimized perturbative mass. Of course, this is only due to fixing from optimization a value of the vari- 
ational mass, and the resulting fixed coupling is similarly a 'variational' parameter. But those variational mass and 
coupling are to be simply replaced in the relevant physical quantity being optimized, giving a relation e.g. between the 
latter and the basic scale A^s of the model, typically given in the MS scheme. Actually, when such RG-improved OPT 



Typically most of our results in refs. [8, 25, 26] depended on further using Pade approximants for the mass dependence, representing 
rather specific ansatzes for the relevant chiral symmetric limit. 
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is applied to the dimensionlcss ratio of the mass gap to A^jg, it is completely cqiiivalcnt to optimizing independently 
with respect to the two mass and coupling parameters. All those properties will become clearer as the procedure will 
be worked out below on a definite model. 

Having thus summarized some previous developments of the ordinary method, the present work is organized as 
follows. Sec. II briefly reviews some known perturbative and non-perturbative results for the 0{2N) Gross-Neveu 
model and sets our conventions, defining also RG quantities to be used later on. In Sec. Ill we present the LDE 
method and the interpolated GN model, supplemented with the new idea on the incorporation of renormalization 
group requirements. For illustrative purposes we first deflne the ingredient of the new RG improved variational method 
on the particularly simple N ^ oo case, where both the exact mass gap is known [9] and its perturbative expansion 
can be derived to arbitrary large order. In this case the new procedure is particularly simple and entirely analytic. 
It displays remarkable convergence properties as the order of the perturbation expansion increases. We discuss it in 
some detail as the generic behaviour of the solutions serves as a very useful guide to the more general case of arbitrary 
N. In section IV we extend this framework to the case of arbitrary N values, examining also as an intermediate step 
the next-to-leading 1/N case. We determine approximations to the mass gap for arbitrary values of N and compare 
those with exact known results to estimate the error of the method. We discuss some difficulties encountered and a 
few variants of the method. Finally conclusions are given in section V. 



II. GN MODEL PERTURBATIVE AND NON-PERTURBATIVE RESULTS 



The 0{2N) Gross-Neveu model [9] is described by the Lagrangian density for a fermion field given by 

£ = § ^) - + i(7^^ + Cc.T. (2.1) 

where \[' is a 0{2N) vector and summation over N is implicit in the above equation. In 1 + 1 space-time dimensions 
the tpk components of ^' for A; = 1, ..A'' represent two-component Majorana spinors. The mass m and coupling 
designate implicitely renormalized quantities and the appropriate counterterm are not shown in Eq. (2.1). 

When m = 0, in addition to the 0{2N) symmetry the theory is invariant under the discrete chiral symmetry (CS) 

*-S-75* , (2.2) 
which is spontaneously broken such that the fermions get a non-zero mass [9]. 

For studying the model Eq. (2.1) in the largc-A^ limit it is convenient to define the four-fermion interaction as 
gQj^N/ir ~ X. Since g^j^ vanishes like 1/A^ we study the theory in the large- A" limit with fixed A [9]. At leading order 
of the 1/A^-expansion, this mass gap is simply 



where in this normalization 



Mjv^oo = Ams (2.3) 



Aj^ = fie (2.4) 



in terms of the renormalized coupling A(/i) in the MS scheme at the renormalization scale fj,. 

Next, the exact expression of the mass gap for arbitrary N >2 values has been calculated [28] from the Thermodynamic 
Bethe Ansatz (TEA). The result is 

^(^)exact ^ (4e)^ ^2.5) 

Ams r[l - 2(jv-l)] 

and also useful is its 1/N expansion at next-to-leading (NLO) order 

M(N) , , l-h21n2-7B , , 

-±J(NL0 1/iV) = l + (2.6) 
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where 7^ is the Euler-Mascheroni constant. 

Independently of those exact resuhs for the massless theory, one can consider in the massive GN 0(2N) model the 
perturbative expansion of the pole mass M in terms of the running mass m = m{p,), which is known at present only 
up to two-loop order: 



1 + 9icr - 7oi) + 9^ {cr + {lo - 7f ^ - crijo + bo)) i + I (70 + 60) i') ] (2.7) 



where we take the normalization: g = g{p) = gQ^^j-K. For the sake of generality we have made explicit in Eq. (2.7) 
the dependence upon the RG beta function i5[g) and anomalous mass dimension ^mig) coefficients, which are known 
up to three- loop order in the MS scheme [29]. (NB the 'MS' indices in Eq. (2.7) indicate that the relevant coefficients 
are scheme dependent, here given in the MS scheme). One defines as usual: 

m = - W - bWg^ ; (2.8) 

and 

7m(5) = 705 + 7r5'+ 71^5'; (2.9) 
and for the 0{2N) GN model one has the particular values in our normalization: 

6o = iV-l; &i = -^; bT = -^{N-l){2N-l); (2.10) 



70 = A^-^; 7r = -|; ir = -YQ^N-\){m-2,). (2.11) 



Also in Eq. (2.7) the finite non-RG perturbative terms are [8, 29, 30]: 



cr=0; cr = {N-\)C-^-^); (2.12) 
The full RG operator reads in this normalization 

which should give zero at some perturbative order g^ (up to terms of higher orders 0{g^'^^)), when applied to a 
physical RG-invariant quantity like the pole mass above. Finally we also give for reference the expression of the basic 
scale Aj^, relevant for instance in Eq. (2.5), and defined in this normalization as[8, 28] 



A„^pe-w^(_-^) ^ (2.14) 

III. THE INTERPOLATED MASSIVE GN MODEL 
A. the LDE method including RG information 

Let us first examine the standard implementation of the LDE proccdinc within the GN model, before supple- 
menting this construction with the new ingredient from RG information. According to the usual LDE interpolation 
prescription [1], from the original four fermion theory, Eq. (2.1), we define the deformed Lagrangian, 

C{m, g,5) = ^{i - m(l - (5)** + 5^{^^f + Cct,s (3.1) 
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so that the new perturbation parameter 5 interpolates between a free massive theory for 5 ~ and the original 
massless interacting theory for (5 = 1. In Eq. (3.1) the connterterm Lagrangian density, Cct.S, has the same 
polynomial form as in the original theory, while the coefficients are allowed to depend on 6. In fact for any given 
ordinary perturbative expression in terms of the original renormalized mass and coupling, in the MS scheme, the LDE 
procedure is implemented simply by the substitutions 

m m(l - S) , g g6 , (3.2) 

into the perturbative series of any relevant physical quantity, and re-expanding the result to order k in the new 
expansion parameter 6. Next, S is set to the value (5=1, such as to recover the original (massless) theory at infinite 
order, while at any finite order k there remains a dependence on m in the LDE of the physical quantity. Considering 
for the latter typically the case of the pole mass, expanded at order k: M^'^\m,S = 1), the standard, mostly used 
prescription is an optimization criterion (principle of minimal sensitivity (PMS) [4]), requiring at each successive 
perturbative orders k: 

■^M^''\m,g,d=l)\m^in = (3.3) 
dm 

thus defining an optimal A;-dependent mass ifi value. Up to now we have described what the standard procedure is 
in most similar studies, up to some variants in the prescriptions to fix the mass^. The main new ingredient that we 
add to this procedure, is to supplement the standard mass-fixing prescription (3.3) with a direct perturbative RG 

information. Namely, at any given order k of the new LDE expansion, we impose, in addition to the mass optimization, 
that the modified LDE series satisfies a standard RG equation at the appropriate perturbative order, i.e. 

RG (m(*=) (m, g,S = 1)^=0 (3.4) 

where the RG operator is defined in Eq. (2.13). Now, combining it with Eq. (3.3) immediately shows that the RG 
equation takes a simpler, reduced form: 

MW(m,5,(5= 1) = (3.5) 

Note also that applying both Eqs. (3.3) and (3.5) completely fixes m = fh and g = g, since one has two constraints 
for two parameters in this case. A further conceptual simplification occurs, when considering directly the required 
ratio M^'^^/Ams-: In fact, since Afj[g{g) satisfies by definition 

Ams^O (3.6) 

consistently at a given perturbative order for P{g), it is easy to show that Eq. (3.3) and (3.5), are completely equivalent 
to the following: 

d /'M''{m,g,6=l)\ d /M'=(to, 5, (5 = 1) 

dm \ K-^ ) ' dg \ Aj^ 

i.e. the procedure is equivalent to optimizing independently with respect to the two parameters of the theory! 




^ For example some other studies use a so-called 'fastly apparent convergence' (FAC) prescription, requiring at order k the order fc -|- 1 to 
vanish. Though this prescription is often simpler (analytically) than optimization with respect to the mass, it needs the knowledge of 
a priori higher orders so is even less practicable than optimization in models where only few perturbative terms are known. Moreover 
the mass optimization is best suited for our additional requirement to incorporate the RG information as we shall see. 
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B. Application to the N oo case 



Let us examine now in some detail how this new RG+OPT prescription works in the leading iV — >■ oo limit of the 
GN model. The N ^ oo limit is particularly suited for analysis because one knows the perturbative expansion of 
the original theory for m 7^ to arbitrary perturbative order. The latter expansion can be cast into the simple form, 
after redefining A = NgQj^/n: 

M{m,X)=m(l + X\n—\ (3.8) 



where m = m{^) and A = are the renormalized mass and coupling in the MS scheme. The simplicity of the 
GN model in the leading large N limit is illustrated in the compact form of Eq. (3.8), containing only the leading 
logarithmic dependence on m, i.e. there are no other next-to-leading logarithmic or non-logarithmic corrections to 
the relation between the mass gap M and the running perturbative mass m. Moreover the correct result M = A.-^ in 
the chiral symmetric limit m ^ can be retrieved directly from algebraic manipulation of Eq. (3.8), noting that the 
latter equation can be rewritten as [8, 27] 

M(m) 



= e^(™/^ns) ^1 for TO , (3.9) 



where we used Ams = fxe^^^^, W{x) is the Lambert implicit function defined as W{x)e^'^^^ = x, r'n = to/A is the scale 
invariant mass, and in the last limit we used the properties W{x) ~ a; for a; ^ 0. (NB clearly the result M/A-^jg = 1 
had been obtained previously by direct calculation in the original massless theory to = [9] , but the previous algebraic 
manipulation makes clear the link with the perturbative RG-resummed form Eq. (3.8) for the massive theory). 
Now assume for a while that we only know at some finite order the pcrturbativcly expanded Eq. (3.8), in terms of 
TO at a given order A*^, and let us examine the results of the modified LDE: after the substitution (3.2) we obtain to 
first order in 6: 

M^^\m,X,S) = m(l-5{l + X\n^)^ . (3.10) 
Now taking 5 — >■ 1, the OPT equation (3.3) gives: 



with immediate solution (assuming A 7^ 0): 



A(l + ln-)=0 (3.11) 



- = e-' . (3.12) 



At this stage, a standard treatment of the LDE+PMS would consist in replacing this optimal mass value within the 
expression of M^^\m, A, 1) and to proceed similarly at successive orders of the expansion. This gives a result for the 
mass with a non-trivial dependence in A very different from the perturbative one. For instance at first order it is 
simply: 

M(i)(m, (y, 1) ^ A TO = A/xe"^ . (3.13) 

However this does not yet give the required relation for M as only function of Aj^. Moreover for finite N, or similarly 
in another more complicated model, it is even less obvious to extract from such result the relevant ratio M/ A^js , due 
to the fact that the OPT equation gives a generally involved gf-coupling dependence, which is not directly related with 
Ams. Let us examine what our new simple prescription gives for the above N ^ 00 series: The (3 function for large 
N is simply given as /3(A) = — A^ (with no higher order corrections) so that the extra RG equation takes the form: 



dX 



M^''\m, X, 6 = 1) = 0, (3.14) 
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which gives at first order when applied on expression (3.10): 



TTi 

mA(l + Aln— ) =0 , (3.15) 



which for m ^ and A ^ gives the non-trivial solution 



A = (ln-^)-i (3.16) 
m 

which reminds of the perturbative expression of the running coupling A for m (the exact running coupling for 

N ^ oo being given by Eq. (3.16) but with m ^ms)- Next combining Eq. (3.16) and Eq. (3.12) gives the result: 

TTI ~ 

In— = -1; A = l, (3.17) 
M 

which upon replacing this solution in the corresponding expression for the pole mass at this perturbative first (5-order, 
Eq. (3.13) gives: 

^"^(-'^1^ = ^) = 1 , (3.18) 

Ams(A) 

i.e. we obtain the exact mass gap result already at the very first order. Eq. (3.17) also implies m = Ams, so that the 
exact running coupling is also obtained. 

One may proceed to successive orders, and in fact at any arbitrary perturbative order 6 we obtain solutions (3.17) 
giving thus the exact mass gap result. This is quite satisfactory, but it is instructive to examine further the behaviour 
of such solutions. Actually, a definite drawback of the optimization prescription is that it involves minimization of a 
polynomial equation of order k in the relevant mass parameter m at perturbative order 5'^, supplemented now by the 
extra RG equation. It is clear that more and more solutions are to be considered when increasing the order and this 
non-uniqueness of the optimized solution may require extra choice criteria. In fact, up to order fc = 3 in the A/' — > oo 
case the (multi)-solutions are degenerate and give uniquely the solution of Eq. (3.17) giving the correct result for the 
mass gap. For instance at order k = 2 one finds 

{1 + XL){1 + L) = (3.19) 

from RG Eq. (3.5), where we defined for convenience L = lnm//U, and 

(1 + A + AL(3 + L)) = (3.20) 

from OPT Eq. (3.3). Although the two equations are now intimately related and are solved together within our 
approach, to make the link with the standard RG properties it is convenient to consider the RG equation as determining 
the coupling A as function of the mass dependence L, while the OPT equation finally determines the optimal mass. 
But note that both possible solutions of the RG equation, the 'standard' RG behaviour A = —L~^, or the other 

sohition L = — 1, lead finally to the same solution of Eq. (3.17) once substitued into the OPT equation (excluding 
the trivial solutions m ~ or/and A = 0). Now things become different starting at third order of the (5-expansion, 
A; = 3, where spurious solutions appear in one or both equations. More precisely the RG equation gives 

(l-|-AL)(l + 2A-(-Ai:(5 + 2L)) = 0. (3.21) 

So, while the previous A = —L~^ is still a solution, having the standard RG behaviour, an extra solution: 

X = -{2 + 5L + 2L'^)-'^ (3.22) 

appears, having clearly the wrong RG-behaviour even for large L. Injecting this very odd solution into the OPT 
equation gives 

(1 + L)3 (7 + 2L) = . (3.23) 
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So even if the expected solution L = ~1 reappears (and recovering again A = 1 from Eq. (3.22)), an extra solution 
remains, that we consider spurious, giving 

L=-l' ^ = -^ (3.24) 

which finally gives M^'^^ (spurious)/ = —354 e"^^^/^'' ~ —1.15 x 10~* which is clearly unphysical. This solution 
may be rejected even if not knowing the right exact result, on the basis that it does not have the correct perturbative 
RG behaviour (moreover giving a negative coupling g). But more generally this illustrates that the higher order and 
non-linear equations implied by the optimization procedure may lead to spurious solutions, as we shall also encounter 
in the less trivial situation of finite N. Also, high order polynomial equations often have unstable solutions, i.e. 
small changes in the coefficients of the higher orders may induce large variations of these solutions. In that case, 
qualitative considerations using the expected RG behaviour plus other criteria may be necessary to remove such 
spurious solutions. 

In summary at the leading order of the l/A''-expansion, we obtain for arbitrary high perturbative order k the 

following properties: 

-the RG equation factorizes to the form 

(1 + AL)/('=)(A,L)=0, (3.25) 

i.e. whatever the form of f^''\X,L), A = —{L)~^ is always a solution. 
-Injecting the latter solution into the OPT equations gives 

(1 + L)*= = , (3.26) 

i.e. giving ^-multiplicity roots of the form Eq. (3.17) thus giving the exact result M^''^ = A-^. 

Moreover, in fact injecting this RG 'perturbative' behaviour solution directly into M^''\m,X, 1) simply gives at any 
fc-order: 

MW(m,A,l)=m (3.27) 

without any extra correction, and having not yet used at this stage the OPT equation determining m. This 
coincidence here between the 'pole' mass Af and the perturbative mass m, meaning there arc no pertiirbative 
corrections, is certainly a peculiar property of the large N limit, but nevertheless a non-trivial result specific from 
this construction, and not obvious from the original form of Eq. (3.8). A similar situation occurs in the large- A'' limit 
of the 0{N) oscillator where one of the solutions of the optimization procedure at each order gives the exact answer 
for the ground state energy with a similar increase in flatness around that optimum [7] . 

-Finally we observe that for any k > 2 the second derivative with respect to the mass is zero at the extrema points, 
i.e. the extrema becomes more and more flat, as intuitively expected. 

To appreciate further the convergence properties of the procedure, we have plotted the surface M*^*^^ (x = A, j/ = In ^) 
for successive values of k in Fig. 1 (using Mathematica [31]). It is striking to sec that this surface becomes flatter 
and flatter around the optimum x = 1, y = —1. This is conflrmed by the expansion of M'*^) around that point: 
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= 1 + 
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- 2)"^ (5a; + 7y + 2) + quintic terms , 




= 1 + 
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- 2)^ (7a; + lly + 4) + order 7 terms 




= 1 + 




-y- 


- 2)^ + order 7 terms . 



At order 2k— 1 and 2k we may conjecture that there is around that point a behaviour of order (fc + 1) in both x and 
y directions, corresponding to the increasingly flat behaviour which the figure suggests. We may see this property in 
the N oo case as an empirical convergence proof of the procedure for any value of A and m in a large neighbourhood 

of A = 1, ln2l = -1. 

All those properties of the large N case will be a useful guide for the less trivial arbitrary case. 

IV. ARBITRARY N CASE 

The modified LDE pertubation given from Eq. (3.2) is now applied on the original perturbativc series given 
by (2.7), only known exactly at two-loop order. In order to examine the eventual improvement and convergence 
properties of the method when increasing the perturbative order, we shall compare the prescriptions at the available 
perturbative orders, namely first and second orders. We will also examine diff'erent prescriptions and approximations, 
also concerning the unknown higher orders to estimate the sensitivity of our results to the latter. 



A. OPT and RG at first order 

After substituting (3.2) into Eq. (2.7), expanding to first (5-order and taking 6 = 1, one obtains 

M^'\m, g, 5 = 1) = -{N - h m g In"^ . (4.1) 

We take the RG Eq. (3.5) at first order for perturbative consistency, i.e. only with the bo dependence entering. At 
this first order, both the OPT and RG Eqs. (3.3), (3.5), or equivalently Eqs. (3.7), respectively take a very simple 
form: 

g{L + l) = (4.2) 

and 

mg{l+gL{N -1))=0 (4.3) 
where L = ln^. The unique non-trivial solution is: 



10 



Note that the sohition of the OPT equation (4.2) is similar to the TV — ^ oo case above, while now the TV-dependence 
enters through the RG solution (4.3). Moreover what is quite remarkable is that just like in the TV — > oc case, the 
RG solution again gives the correct RG behaviour of the running coupling for /n ^ m: g{ij) = (6o In ^)^^, since Eq. 
(4.4) also implies m = Ai^, noting that at first RG order, = p,e~^/^^^~^^^\ Now, the correct comparison with 
the exact result Eq. (2.5) implies to use the very same normalization for entering Eq. (2.5), defined in (2.14), 
and to substitute for g the optimal solution (4.4): we thus obtain for the ratio of (4.1) and A^s- 

- ' 2(1 - ^771^^^) ■ (4.5) 



Ams N-1 y 2 (AT - 1) 

This gives 3/2, ~ 1.38, ... for A/' = 2, 3, ... to be compared with the exact result Eq. (2.5), e.g. ~ 1.8604, 1.4819, . . . 

for A' = 2, 3, Apart for A' = 2 for which it is a rather poor approximation of the exact result, this is already quite 

reasonable for A^ = 3, and better for larger A" since the discrepancy decreases as 1/A^. (Note also that Eq. (4.5) is 
zero if analytically continued to A^ = 3/2, consistently with the exact result (2.5)). More interestingly, the OPT first 
order result compares even much better with the next-to-leading order (NLO) in 1/A' expansion of the mass gap: 
expanding to first 1/A' order Eq. (4.5) one finds: 

M(i), , l + ln2 , , 

^(NL0 1/Ar) = l + ^, (4.6) 

to be compared with Eq. (2.6): the difference, (ln2 — 7b)/(2A') ^ 0.058/A' gives a relative error of Eq. (4.6) with 

respect to (2.6) of only ^ 2% for = 2 and much less for N > 2. This result at NLO in 1/A" can be attributed to the 
fact that the 6o first order coefficient in the RG equation is actually the only contribution, all other RG dependence 
being 0{1/N^), and thus the first (5-order RG solution (4.4) turns out to give the correct complete RG dependence 
at this 1/N order. 



B. (5-expansion at second order 



We consider now the second order in the (5-expansion applied to expression (2.7) after substitution (3.2). At this 
order, the RG and OPT equations can still be managed analytically. Solving consistently at order 6"^ we obtain 
respectively for the OPT Eq. (3.3): 

g [48 + 5 (-21 + 48Af + 12L(-7 - 3L + 4(3 + L)N) + iir^) ] = , (4.7) 

and for the RG equation (3.5): 

m c/2 [12 (5 + 6L - 8(1 + L)N) - 2g{N - 1) (-21 + 12L(-1 - 3L + 4(1 + L)N) + A-k'^) (4.8) 
+g'^{N - 1) (-9 + 12L(-1 - 3L + 4(1 + L)N) + A-k'^) ] = . 

(which again also exhibit trivial solutions g = Q and m = that we of course ignore). 

Let us consider first the case A" = 3 as a typical moderate value {i.e. sufficently far from large N). Solving the two 
coupled RG and OPT equations (4.7), (4.8), one finds five different solutions: three give in fact complex-valued m, g 
and M^'^\m,g), that we reject at the moment as a priori unwanted 'spurious' solutions, and the two remaining real 
solutions are 

g ~ 1.69, L ~ -2.52 ^^^i^'^^ ~ 3.65 (4.9) 
g ~ 0.41, L ~ -1.74 ^ 1.503 

so that the second solution is very close, about 1.4% above the exact result ~ 1.48185. The other real solution is far 
away, and we now argue that it can be considered spurious. First, the relatively larger absolute values of g and L for 
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TABLE I: Combined OPT+RG results at second (5-order, selecting solution (4.10), for the GN mass gap for different AT, with 
the corresponding values of optimal mass parameter rh and coupling g, as compared with Af^^^'^'/A^jg-. 



N 


M(^' (m,<j) 


L = In 22 


rh 


9 




% error 




MK 




2 


1 69 ± 09 i 


—2.28 ± 0.18't 




1.08 ± 0.34i 


1 86038 




3 


1.503 


-1.74 


0.70 


0.41 


1.48185 


1.4 % 


4 


1.338 


-1.57 


0.81 


0.27 


1.3186 


1.5% 


5 


1.252 


-1.47 


0.89 


0.20 


1.23668 


1.2% 


8 


1.143 


-1.33 


0.92 


0.12 


1.1330 


0.88% 


10 


1.110 


-1.28 


0.95 


0.094 


1.10285 


0.65% 


100 


1.0099 


-1.019 


1.005 


0.0099 


1.00916 


0.07% 



this solution may indicate that it conies from an excessive influence of the higher powers in the polynomial equations, 
which powers are not to be trusted at this order. Furthermore, a closer inspection of Eqs. (4. 7), (4. 8) is easy to do 
for any A'', since the RG equation is second order in g: expressing the two resulting solutions for g{L), one reahzes 
that only one solution has the correct 'perturbative' RG and large N behaviour: expanding the solutions for large \L\ 
gives: 

where the first leading term exhibits tlie correct large A'^ and leading logarithm expression for the running coupling: 
~ (^oln^)~^ for /X ^ m (provided that m/Ans remains of 0{1), similarly to the leading and next-to-leading 
1/A'' orders above), while the other solution gives 

1 2 -I- A/' — AN'^ 1 

'''"'^ - ' + (¥31)1 + 2(iV -\f (4iV - Z)L' + "'l?' C^") 

in clear conflict with both the known large A^ behaviour 5 ~ 1 /N and a standard perturbative RG behaviour of tlic! 
coupling. According to the lesson of the large A^ limit, we are lead to select uniquely the first solution, which indeed 
is the one giving M/ kyis{N = 3) ~ 1.503. The same is observed for any A^ > 3 values: according to the RG behaviour 
criteria, we uniquely select the first solution and this gives the result collected in Table I. We also give for illustration 
the corresponding values of fh/h.-^: one can remark that for arbitrary N >3, the solutions are such that all m/Ay[s 
values remains of order 1 and regularly tends towards 1 for increasing A^, which is an a posterori crosscheck of the 
right perturbative RG behaviour of the solution (4.10). Moreover, the absolute value of the second derivative with 
respect to the mass is smaller for this solution, which is an additional indication of the expected behaviour similarly 
with the large N above results. Also, in analogy with the A' ^ 00 case, we can study for arbitrary N the surfaces at 
second (5-order, M^'^\x = g,y = lnm//i): let us just mention that from a direct minimization in the {x,y) plane the 
same values for A^ > 3 of Af^^^ /A^^ as those in Table I are recovered, and that these surfaces exhibit an increasingly 
flat behaviour around those extrema for increasing N. 

A rather embarassing problem, however, occurs with lower values of N: actually, for AT = 2 we find only complex 
solutions when solving Eqs. (4.7,4.8). In table I we put the result with the smallest imaginary part, but even when 
taking the real part only it is not very close to the exact result, as compared with A^ > 3 values obtained by the 
same procedure. We can examine this problem more closely by zooming on A^ values: since the TBA results arc 
formally defined for any real N > 1, we can always continue this expression for non-integer A^ values and compare 
to our approximation procedure, which is evidently also well-defined for any real N. The results are given in Table 
II for a few representative 2 < A' < 3 values. Our second-order approximation remains excellent for values N < 3, 
until the point when solutions becomes complex, which occurs at A" ~ 2.1 approximately. This refiects the fact 
that the optimization prescription Eq. (3.3) alone already involves a polynomial equation of order k for Inm, and 
the coupled OPT-RG equations (3.7) become non-linear. As a result more and more solutions, some of them being 
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TABLE II: Same as in Table I for some representative real 2 < iV < 3 values. 



N 




L = 111 ^ 




-1/, , „, t 


Yi error 






2.1 


1.72±0.008i 


-2.26 ± O.OSz 


1.08±0.15i 


1.816 




2.2 


1.800 


-2.37 


1.31 


1.768 


1.8% 


2.4 


1.679 


-1.96 


0.642 


1.677 


0.1% 


2.6 


1.614 


-1.87 


0.537 


1.600 


0.9% 


2.8 


1.554 


-1.80 


0.466 


1.535 


1.2% 



eventually complex, are to be considered when increasing the order, and it is not much suprising that at second 
order the occurence of complex solutions also depends on the value of A''. In a previous work [19] we had proposed a 
rather simple way out for this problem with a generalization of the PMS criterion as performed on the LDE series, 
which turns out to lead to a drastic reduction of physically acceptable real optimization solutions at each successive 
perturbative order. The modification is to introduce extra variational parameters within the interpolating Lagrangian, 
starting here at the relevant order two with one more parameter a introduced as 

m^mil-S)" . (4.12) 

Next, the OPT criterion (3.3) is generalized by requiring both the standard optimization (3.3) and 

52m(2)/9to2=0 (4.13) 

giving a system of two equations to be solved simultaneously for a and m, such that one can make real a pair of 
solutions which originally had a relatively small imaginary part for the standard interpolation with a = 1. At higher 
orders, the generalization is easily done with additional parameters and additional vanishing of higher derivatives 
of M(m). Such a modified prescription has been applied successfully e.g. to the calculation of the shift in critical 

temperature of the Bosc-Einstcin condensate (BEG) (where the knowledge of high perturbative orders leads to a 
similar problem of complex OPT solutions), with excellent agreement with lattice results [19]. 

There are of course other possible perturbatively equivalent ways to introduce such an extra variational parameter, as 
long as the only constraint is that the modified Lagrangian still interpolates between the free field (massive) theory 
for 5 = and the original (massless) theory for S = 1. The non- linear (5-expansion with a simple exponent form 
in Eq. (4.12) gives algebraically simpler subsequent optimisation and RG equations (3.7), and is inspired from a 
prescription [19] suited to the D = 3 0{N) g^'^ model case, where it is closely connected to other prescriptions, intro- 
ducing explicitly the critical exponent related to the anomalous mass dimension within a power-modified interpolating 
mass [2, 16, 18], which was argued to drastically improve the OPT convergence. (NB indeed the exact value of this 
critical exponent for the D = 3 0{N) model at the next-to-leading 1/N order had been obtained empirically in 
[19] by solving the analog of Eq. (4.13) together with the analog of the OPT Eq. (3.3). In the present D = 2 GN case 
however, this modified OPT prescription is to be viewed as an essentially algebraic trick to force one of the solutions 
to be real for A'' ~ 2, with a priori no particular physical meaning of the extra variational parameter a). 

Using Eq. (4.12) instead of (3.2), expanding (2.7) at order 5"^, and applying Eqs. (3.7) and (4.13) to solve for 
o, m and g gives the result in Table III for some representative values of 2 < A'' < 3. We recover in this way 
real results for N ^ 2 which are very good for N < 2.5. However, when N is increasing, results are somewhat 
worse than the original ones in Table II from the simplest OPT prescription. Of course, as long as we directly 
obtain from the simplest prescription real solutions for any A'' > 2.1 values, it is not needed to appeal to the more 
elaborate interpolation with an extra parameter a: thus we obtain in this way an overall controllable prescription 
for any N. Moreover even if we would not know the exact solution, it is easily seen in the present case that the 
solutions for a 7^ 1 progressively depart from the expected 1/A^ and perturbative RG behaviour for increasing N values. 
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TABLE III: Combined extended OPT+RG results at second (5-order for the GN mass gap for some representative real 2 < iV < 3 
values, with corresponding values of the optimal parameters a, m, g, and compared with M'^^"'''^ / h.-^. 



N 


M(^) (m,ij) 


a 




9 




% error 


^dj^ 


^iij^ 


2 


1.907 


1.986 


-3.815 


0.187 


1.86038 


-f2.5% 


2.1 


1.812 


1.931 


-3.716 


0.170 


1.816 


-0.2% 


2.2 


1.734 


1.886 


-3.637 


0.155 


1.768 


-2% 


2.4 


1.614 


1.818 


-3.516 


0.133 


1.677 


-3.9% 


2.6 


1.526 


1.77 


-3.428 


0.116 


1.600 


-4.8% 


2.8 


1.458 


1.734 


-3.36 


0.103 


1.535 


-5.3% 



C. Approximate schemes at second order 

In the previous second J-order calculation we required the exact two-loop order RG equation (4.8) to hold, which 
results in a second order equation in g (after eliminating the trivial solution g = Q), thus responsible for the occurence 
of two solutions in Eqs. (4. 10), (4. 11), one of those being spurious according to the RG behaviour, or complex for small 
N values. In fact one can think of two alternative ways to avoid this problem from the beginning, while still using 
the RG perturbative information, noting that the purely perturbative information at this available second order does 
not require the RG equation to hold exactly. Firstly, standard RG invariance only requires that the full RG operator 
(2.13) when applied to (2.7), gives a remnant term of 0{g^), rather than being exactly zero, since those higher order 
terms would be compensated by higher order (3-loop) terms in the pole mass itself, which are not available. We 
can impose a similar requirement in the case of the ^-expanded pole mass expression, truncating perturbatively the 
reduced RG Eq. (3.5) instead of requiring it to hold exactly. Secondly, quite similarly, the expansion of the exact 
second order solution of g{L) for large \L\ in Eq. (4.10) formally gives an infinite series for the leading, next-to-leading, 
etc. .. logarithm L dependence (LL, NLL). But, the genuine NNLL etc. ..orders also need the knowledge of higher 
(3-loop and beyond) RG behaviour, so that in ignorance of the latter, wc may truncate Eq. (4.10) at the NLL Level. 
Those two possible approximation schemes for the RG dependence are not numerically equivalent, as we shall see, 
and such variants may also provide indirectly an estimate of a theoretical error of the method at this second order. 

1. Perturbative truncation of the RG equation 

In analogy with the ordinary RG properties, we only require a third-order tnmcated reduced RG Eq. (3.5) (or 
equivalently Eq. (3.7) for the RG part), when applied to the variational (5-expanded pole mass, since the fourth order 
actually contains RG terms that would cancel with three-loop terms in the pole mass. This third order truncation of 
Eq. (4.8) nevertheless keeps a dependence on the two-loop RG coefiicients bi and 71, but now gives a simple linear 
equation (omitting of course the trivial solution g = 0) thus with a unique solution: 

^~N+{^-N)L-g [21 - An^ + 12L (1 + 3L - 4(1 + L)N)] = (4.14) 

which moreover is easily checked to have the correct large and RG perturbative behaviour: g ^ — l/((iV — 1)L) 
for large \L\. Combining this solution with the OPT equation gives a third order equation for i, with a unique real 
solution, which moreover is such that fh/Aj^ is very close to 1 for any N > 2 values. This gives the results for 
arbitrary N given in Table IV. We note that those are the closest results to the exact mass gap, with the relative 
error being less than 1% for any N >2, and often well below the percent level. 
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TABLE IV: Combined OPT+ third-order truncated RG at second (5-order for the GN mass gap for different N, with the 
corresponding values of optimal mass parameter m and coupling g, as compared with Af^'^^'/A^jg-. 



iV 


M<2>{m,9) 


L = In 22 


m 


g 




% error 


A 


A 




o 

Z 


1 Q A A '7Q 

i.o44/o 


1 KAQ 

-i.o4o 


n CQ /I 


n QQ /I 


i.oDUoo 


-U.oD/o 


3 


1.48460 


-1.313 


0.958 


0.443 


1.48185 


-h0.19 % 


4 


1.32627 


-1.209 


0.978 


0.307 


1.3186 


-1-0.58% 


5 


1.24448 


-1.138 


0.984 


0.239 


1.23668 


-hO.63 % 


8 


1.13956 


-0.962 


0.979 


0.158 


1.133 


-F0.58 % 


10 


1.10896 


-0.952 


0.981 


0.123 


1.10285 


-hO.55 % 


100 


1.00993 


-0.958 


0.998 


0.0106 


1.00916 


-1-0.077 % 



TABLE V: Combined OPT-|- g{L) dependence truncated at NLL level, at second J-order for the GN mass gap for different N , 
with the corresponding values of optimal mass parameter fh and coupling as compared with M"^'^'^* / A^. 



N 




M 


m 






% error 






9 




2 


1.8663 


-1.77 


0.69 


0.82 


1.86038 


+0.3% 


3 


1.503 


-1.55 


0.85 


0.41 


1.48185 


+ 1.4 % 


4 


1.337 


-1.45 


0.92 


0.27 


1.3186 


+ 1.4% 


5 


1.252 


-1.39 


0.90 


0.21 


1.23668 


+ 1.2 % 


8 


1.143 


-1.29 


0.96 


0.12 


1.133 


+0.9 % 


10 


1.1106 


-1.26 


0.95 


0.095 


1.10285 


+0.7 % 


100 


1.0099 


-1.07 


0.997 


0.0095 


1.00916 


+0.07% 



2. Truncation at the next-to-leading logarithm order approximation 

Next wo examine here the results of a different kind of truncation, as discussed above, namely truncating the 
right behaviour solution Eq (4.10) for large L at the NLL level, invoking that in a standard perturbative framework, 
higher NNLL etc. . . orders need higher (three-loop) order information that we do not take into account in Eq. (2.7). 
Combining this with the OPT Eq. (4.7) gives the result in Table V. Those results are reasonably good, at the percent 
or so level, but not as good as the previous ones obtained from truncating directly the RG equation. 

D. Higher order estimates and stability 

Our results above give a clear evidence of numerically fast convergence of this RG-improved implementation of the 
OPT method. In the absence of a rigorous convergence proof, we can still try to examine higher order estimates and 
the sensitivity and stability properties of our results with respect to higher orders. The known three-loop RG beta 
function and anomalous mass dimension in Eqs. (2.10,(2.11) may be used, giving all the logarithmic dependence at 
order in an extension of Eq. (2.7): 

Mf-i"""^ = m [1 + 5(ci - 7oi) + 9^ {cW + {ll - 7p - ci(7o + 6o)) i + ? (lo + &o) i') (4.15) 

+5^ (aso + asi + x + C3) ] 

where the general expressions of the LL, NLL, and NNLL terms can be derived by applying RG properties: 

aso = -70/6(70 + + 26o) (4-16) 



031 = ci(7o + 6o)(7o + 26o)/2 + 71(70 + &o) + lo{bi - 3/2&o7o - 7o) 



(4.17) 
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TABLE VI: Third order estimates for combined OPT+ RG for the GN mass gap for different N, as compared with M"""''* / h.-^, 
with corresponding values of the unknown three- loop non-RG coefficient C3. 



TV 


M(^) (fh,g) 


M<-^>(m,9) 


M(^) (m,ij) 




(C3 = -bl) 


(C3 = 0) 


(C3 = bi) 




2 


1.845 


1.898 


1.96±0.03i 


1.86038 


3 


1.411 


1.4666 ±0.0005i 


1.54 ±0.04i 


1.48185 


5 


1.192 


1.2388 ± 0.008i 


1.27±0.04i 


1.23668 


100 


1.006 


1.0108 iO.OOli 


1.014 ±0.004i 


1.00916 



(132 = -02(70 + 26o) - ci(6i + 71 - 70(70 + bo)) - 7o + 27071 - 72 (4-18) 

with in the GN model ci = and C2 and the bi and 7^ coefficients are given above. But the three-loop non-RG 
coefficient C3 in Eq. (4.15) is presently unknown. To have a very conservative estimate of the latter, we do not even 
assume it to be positive, and vary it between < \c^\ < \c3\max where for |c3|maj, we take a "theoretically-inspired" 
maximal value, namely the coefficient given by the leading renormalon [32] factorial behaviour at large perturbative 
order. This factorially growing behaviour of the purely perturbative GN pole mass has been estimated from graphs 
at the next- to- leading 1/N order in ref. [33], and is very similar to the well-known leading renormalon of the QCD 
pole mass [32]. In the present normalization, it gives for the coefficient at perturbative order k: 

|cfc|~l(fc-l)!6S-^ (4.19) 

Taking three representative extreme values, C3 = and from Eq. (4.19) C3 = ±6q = ±(A'' — 1)^ as a crude estimate 
for the 3-loop coefficient C3 in Eq. (4.15), we applied the OPT and RG equations on the third order ^-expansion of 
Eq. (4.15). This gives an OPT equation cubic in L and a RG equation quartic in g, thus as expected the combined 
higher-order equations (3.7) have many solutions, most of those being complex and some solutions being quite 
unstable. But similary to the previous case at order 2, it is not difficult to select the much fewer solutions having the 
correct perturbative and large- behaviour, with the results given in Table VI. For those well-behaved solutions, the 
complex results have reasonably small imaginary parts, but in this case we did not attempt to refine this analysis to 
recover real solutions with the above alternative prescriptions, since it is anyway based on a crude estimate of an 
actually unknown coefficient. 

Overall these results can be considered very stable: a large variation of the third-order coefficient, — 6q < C3 < 6g, 
leads to very reasonable changes in the final result. Moreover even for the extreme lower and upper bounds ±|c3|TOaa;) 
the results bound very well (and remain rather close to) the exact mass gap. Note that the renormalon-inspired 
coefficient IcsItoox grows fastly with N, since 60 = 0{N), but this does not prevent the mass gap result to be 
close to the exact one for sufficiently large N. While some among the multi solutions of Eqs. (3.7) are clearly 
unphysical or can reflect instabilities, the reason is that selecting the OPT -|-RG solution giving the correct large 
N behaviour for the coupling, namely g ~ 1/-/V, implies that the g^ ~ 1/A''^ compensates for the renormalon 
jcsl ~ N'^ behaviour for large N. We also checked finally that there are intermediate values of [csj < b^ such that 
a well-behaved solution is real and very close to the exact mass gap (for instance this happens for C3 ~ 0.61 for A'' = 2). 



V. CONCLUSIONS AND OUTLOOK 

We have examined a very simple additional ingredient to the standard variational OPT approach, in order to 
supplement it with RG perturbative information. The procedure is very transparent and largely analytical at second 
order. It has as immediate consequence to fix both the mass and the coupling such that there are no free parameters 
for the massive GN model typically. Those variational mass and coupling are to be simply used to evaluate the 



16 



physical quantity to be optimized, giving a physically relevant relation e.g. between the latter and the basic scale 
A—. Moreover it is completely equivalent to optimizing independently simply with respect to the mass and coupling 
parameters. 

Although our results are numerical and not a proof of convergence of the method, by adding the RG behaviour 
content of the theory in this most simple way we obtain successful approximations of the exact mass gap below or at 
most at the percent level for arbitrary N values, using only the purely perturbative information. 

One may wonder why the GN model behaves so well, in comparison e.g. with the BEC model where the convergence 
appears slower [14, 18, 19], despite its super-renormalizable properties implying that only the mass has non-trivial 
RG properties. Of course the GN model for large is a particularly simple theory, where the original pcrtm-bative 
coefficients do not exhibit the standard factorial divergences of most field theory perturbative series. But for arbitrary 
N the perturbative relation between the pole and Lagrangian mass, as well as RG properties, are very similar to those 
in more complicated theories, such as QCD typically. In fact, even the large-order behaviour of the coefficients of this 
perturbative series in the massive GN case have a similar, badly factorially divergent, (renormalon) behaviour [33]. In 
that respect the results obtained here with a fast numerical convergence at the percent level or less, and controllable 
from understood RG behaviour, are to be considered very encouraging for similar analysis in more involved D = A 
theories where D = 2 exact S-Matrix and TBA results are inapplicable, and where the mass gap and other related 
non-perturbative quantities are not known and usually not expected to be derivable from first principles. This will 
be the object of future works. 
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FIG. 1: The surfaces M^-^^ as a function of x = A and y = In for A: = 1 , 2 , 3 , 4 (from top left to bottom right) 



